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SEACON  II  STRUMMING  PREDICTIONS 


I . INTRODUCTION 


Seacon  II  (1)  is  an  ocean  based  cable  array  being  used  to  validate 
various  array  motion  models  (2,3,4).  The  array  essentially  consists  of 
three  riser  cables  and  a horizontal  delta  as  schematically  shown  in 
Figure  1.  The  cables  comprising  the  delta  see  uniform  currents  over 
their  lengths  and  hence  arc  subject  to  strumming.  As  is  well  known 
(5,6)  strumming  leads  to  a virtual  increase  of  the  steady  state  drag 
coefficient  of  a cable  over  that  value  measured  for  a stationary  cable. 
Since  this  coefficient  is  a basic  parameter  in  all  array  motion  models, 
an  accurate  knowledge  of  its  value  is  required  to  validate  the  models. 


In  this  report,  the  strumming  behavior  of  the  cables  comprising  the 
Seacon  II  delta  is  examined.  Using  recently  developed  strumming  model- 
ing techniques  (6,7,8)  the  amplitudes  and  frequencies  of  strumming  are 
predicted.  These  amplitudes  and  frequencies  are,  in  turn,  used  to  pre- 
dict the  virtual  increases  of  the  steady  state  drag  coefficients  of  the 
delta  cables. 


II.  EQUATIONS  OF  MOTION  FOR  A TAUT,  NON-UNIFORM  CABLE  WITH  ATTACHED 
MASSES 


We  wish  to  develop  the  equations  of  motion  for  a taut,  non-uniform 
cable  with  attached  masses.  The  situation  is  as  illustrated  in  Figure 
2.  In  this  figure,  s is  the  measure  of  length  along  the  cable  ranging 
from  zero  to  the  total  length  of  cable  L.  The  mass  of  the  cable  per 
unit  length  is  denoted  by  m(s)  and  the  displacement  of  the  cable  from 
its  horizontal  equilibrium  position  is  given  as  y(s,t).  The  concen- 
trated masses  are  given  as  M^,  M^,  ...  and  their  positions  along  the 

cable  represented  by  s^,  The  equilibrium  tension  in  the 

cable  is  given  by  T which  is  considered  large  enough  to  1)  hold  the 
cable  essentially  straight  in  equilibrium  and  2)  be  unaffected  by  small 
cable  motions.  The  above  assumptions  on  T are  equivalent  to  the  taut 
cable  approximation. 


Hence,  the  motion  of  each  cable  segment  between  the  concentrated 
masses  is  defined  by  the  taut  string  equation 


where  t denotes  time.  To  obtain  the  conditions  at  the  concentrated 
masses,  we  must  recall  that  the  cable  has  no  bending  rigidity  and. 
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discontinuities  in  slope  may  occur.  This  situation  is  illustrated  in 
Figure  3 where  F represents  the  force  due  to  the  concentrated  mass.  On 
equating  the  sum  of  the  forces  at  the  mas3  point  to  zero,  one  condition 
at  the  concentrated  mass  is  obtained  as 


T r?  I - 1^  | ~]  + F = 0 

1 os  os  ' J 


The  second  condition  at  the  concentrated  mass  is  continuity  of  deflec- 
tion. Since,  by  D'Alembert's  Principle,  the  force  F due  to  the  mass  is 
given  by 

F = -M  d2y/*t2  (3a) 


equation  (2)  can  be  rewritten  as 


I 

L*s  ' 


- I2  I - M S2y/*t2  * 0 

os  J 


III.  NORMAL  MODES  AND  NATURAL  FREQUENCIES 


We  proceed  now  to  obtain  expressions  for  the  normal  modes  and 
natural  frequencies  of  a taut  cable  with  attached  masses  for  the  case 
in  which  the  mass  per  unit  length  along  the  cable  is  piecewise  constant. 
For  this  case,  it  is  notationally  convenient  to  treat  each  point  along 
the  cable  at  which  the  mass  per  unit  length  changes,  but  which  is  not 
coincident  with  a concentrated  mass,  as  a point  at  which  a concentrated 
mass  having  M = 0 is  located.  From  equation  (3b),  we  see  that  at  such 
a point  the  cable  slope  remains  continuous  which  is  the  correct  transi- 
tion condition. 


The  notation  indicated  in  Figure  4 is  adopted.  To  obtain  the 
normal  modes  and  natural  frequencies  of  the  cable,  solutions  to  equation 
(1)  are  sought  as 

y = tCs)  e (4) 

where  i|r  (s)  is  the  normal  mode  corresponding  to  an  eigenfrequency  tu. 

From  equation  (1)  and  the  adoped  notation,  i|r(s)  has  the  form 


Pi  / mi 

= A^  sin  y — u)  s + Bi  cos  — ui  s for  s^ 


« s ^ 


Applying  the  conditions  of  deflection  continuity  at  the  concentrated 
mass  points,  plus  the  zero  deflection  conditions  at  the  cable  end-points, 
results  in  relations  between  the  coefficients  given  by 


A^  sin 


PT  H\ 

M x 9i  + Bi  cos  1 st 


Ai+1  sin  1 — oo  + B1+1  cos  1 — — » s.  , i = 1 I (6b) 


AI+1  Sin 


V to  L + Bt  . , cos  ^ 


to  L = 0 


The  force  balance  conditions  of  equation  (3b)  yield  additional  relations 
between  the  coefficients  given  by 


V ^ fA.+1  cos  CO  Si  - Bi+1  sin  V~^±i  to  Si  j 

J mi  r * J mi  _ . J mi  1 

-1  — A.  cos  1 — to  s.  - B.  sin  l — t«  s. 

Tl  i T l l ? T i J 


~ fo  |~A  sin  V 
r I l+l  T 


— to  Si  + Bi+1  cos 


m 


r“sij  = 0> 


i = 1,  ....  I (7) 

Equations  (6b),  (6c),  and  (7)  yield  21+1  equations  for  the  21+1  coeffi- 
cients Ai,  ^2*  ®2*  •••»  Ax+1*  ^1+1  * 

These  equations  can  be  written  in  matrix  form  as 

Q C = 0 (8) 

where  C is  a 21+1  column  matrix  having  the  elements 

C - A^,  Aj , Bj » A^»  Bj,  • • • , ^i+x’  ®x+l 

and  q is  a square  matrix  having  rank  2X+1.  The  non  zero  elements  of 
Q are  defined  by 


J mx 

Ql,l  = " Sin  ’ T~  00  S1 

J “ 1 J ml 

q2,i  s^rco8T“,si 

n i J mi+l 

«2i-l,  2i  = sin  ^ — to  Sl 





M 

Jm  i+1 

!2i,2i  T ^ 8in 

1 T 

Vm  1+1 

2i+l,2i  ' 

__  U 

T 

’ Si+1 

- . V mi+i 

cos  " 

mi+l 

!2i+2,2i  T 

T 

Jm  i+1 

■2i-l,2i+l  COS 

— — U) 

T 

Si 

M. 

■2i,2i+l  = r ^ C0S 

% i+1 

’ T 

Jn- i+1 

Si+1 

2i+l,2i+l  1 

m W 

T 

sin  ^ 

HvT 

2i+2,2i+l  ’ T 

T 

i+1 


“ Si+1 


(9e) 


(9f) 


(9g) 


“>  si  (9h) 


(9i) 


In  equation  (9c-9j),  i = 1,2  ...  , I,  The  calculated  elements  Q 
and  Q2I+2  2l+l  are>  however>  ignored  in  forming  Q . 


( 9 j ) 

21+2,21 


The  eigenfrequencies  u)n,  n = 1,  2,  ...  of  the  cable  are  determined 
as  the  solutions  of 


det  ^ Q (m)  = 0 


(10) 


I 


' 


The  corresponding  elements  of  the  eigenvector  C,  (u)  ) are  then  found 


from 


1 _ / 

i , n n 


"l,n 


K-l,n 
"K+ 1 , n 


'21+1, n 


= -c, 


K,n 


Q1,K  (u)n} 


QJ-l,K(uJn) 

QJ+l,K(u)n) 


Q2I+l,K(u’n 


(ID 


where  i?  is  a square  matrix  of  rank  21  obtained  by  eliminating  the 
J row  ana’K  column  from  Q . The  value  chosen  for  C is  arbitrary. 


4 


j i 


i 


i 


l 


IV.  CHARACTERISTICS  OF  THE  SEACON  II  DELTA  CABLES 


The  physical  characteristics  of  the  three  Seacon  II  delta  cables 
are  tabulated  in  Table  1 (9).  In  this  table,  d is  the  diameter  of  the 
cable.  The  tension  in  the  cable  is  taken  as  the  average  of  the  maximum 
and  minimum  tensions  along  the  cable  in  the  no-current  configuration  of 
the  array.  Typically,  the  tensions  along  a cable  varied  less  than  10% 
so  that  the  taut  cable  approximation  is  valid. 


Using  the  information  in  Table  1 in  conjunction  with  equations  (9), 
(10),  and  (11),  the  normal  modes  and  natural  frequencies,  up  to  a maxi- 
mum frequency  of  25  rad/s,  for  each  of  the  three  Seacon  II  delta  cables 
were  calculated  and  stored  on  tape.  Also  calculated  and  stored  on  tape 
for  each  normal  mode  and  natural  frequency  was  the  modal  scaling  factor 

I defined  by 
n 


I = 
n 


'^0 


i|f2  ds 
n 


(12) 


V.  STRUMMING  FREQUENCY  AND  AMPLITUDE  CALCULATIONS 

The  relative  orientation  of  the  Seacon  II  delta,  referred  to  an 
x-y  base  system,  is  shown  in  Figure  5 (9).  In  this  figure,  V denotes 
the  current  magnitude  at  the  height  of  the  delta  above  the  bottom  and 
9 denotes  the  current  direction,  referred  to  the  x-y  base  system,  at 
tKis  height.  Since,  in  actuality,  the  delta  is  slightly  skewed  to  the 
horizontal  and  also  possesses  some  sag,  an  average  height  above  bottom 
of  the  delta  of  2460  ft  is  assumed. 


The  naturally  occurring  shedding  frequency  u>  from  a delta  cable 
is  given  by  the  Strouhal  relation  employing  the  normal  free  stream 
velocity  component  as 

tv  = 2nS  (V/d)  | sin  («-•»)  I (13) 

s v 

Here,  S =^0.21  is  the  Strouhal  number.  In  general,  the  cable  responds 

in  its  n mode  when  uu  ^ uu  £ 1.4  id  with  the  maximum  response  occuring 

for  a)  1.2  u)  . For  the  lower  modes  of  cable  vibration  this  modal 
s n 

response  is  pure;  however,  for  the  higher  modes  of  vibration,  severa^ 
modes  participate  in  the  response  (10)  since  the  bandwidth  for  the  n 
modal  response  overlaps  adjacent  modes.  Little  is  known  about  this 
latter  type  behavior,  and,  in  the  spirit  of  engineering  approximation, 
we  assume  pure  modal  responses  for  the  Seacon  cables. 


To  choose  the  mode  n of  the  cable  response  for  a given  Strouhal 
frequency,  uug,  we  use  the  equation 


U)  5 (D  < 0) 


n+1 


(14) 


The  maximum  amplitude  distribution  Yj^ts)  along  the  cable  in  this 
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response  mode  is  given  by  (6,7) 

Ws)/d  ' W ln  I V“>  I (l5> 

and  the  frequency  of  the  response  is  approximately  u>  . In  equation 
(15),  ^ is  a dimensionless  modal  response  amplitude  and  I is  de- 
fined byequation  (12).  n 

As  has  been  shown  in  refs.  (6,7,8),  the  value  of  A^  is  uniquely 
determined  by  the  value  of  a response  or  stability  parameter  S which 
is  itself  related  to  certain  physical  properties  of  the  strumming  struc- 
ture. The  relationship  between  and  is  shown  in  Figure  6.  The 

legend  for  the  data  points  can  be  Tound  in  ref  (7).  The  solid  line  in 
this  figure  is  given  by  the  formula 


W U29/  (l  + 0-43  v3': 


which  is  a least  squares  fit  to  the  experimental  data  points. 

The  response  parameter  SG  for  a structure  strumming  in  its  n1 
mode  is  defined  by 


S„  = 2nSf!  k 

G S 

where  the  reduced  damping  kg  is  given  by 

k = TTc/2pvp 
s n 


(17a) 


(17b) 


Here,  v is  the  kinematic  viscosity  of  the  fluid,  p is  its  density,  B 
is  the  vibratory  Reynolds  number  given  by 

B = w d2/4v  (18) 

n n 

2 

and  c (lb  s/ft  ) is  the  viscous  damping  constant  per  unit  length  of  the 
structure  at  frequency  u)  as  measured  in  the  still  fluid.  This  damping 
constant  thus  includes  both  structural  and  fluid  effects.  Recently, 
Skop,  Ramberg,  and  Ferer  (11)  have  shown  that  the  fluid  contribution 
to  c,  Cj.,  is  given  by 

cf  = 4.5  rrpv  bJ  (19) 

and  Ramberg  and  Griffin  (12)  have  found  that,  for  the  Seacon  type 
cables,  the  structural  contributions  to  c are  negligible  compared  to 
c ^ • Substituting  equation  (19)  into  equation  (17b),  kg  for  the  Seacon 
delta  cables  is  found  to  be  uniquely  defined  in  terms  of  the  vibratory 
Reynolds  number  of  the  response  as 


k = 4.5  rr2/2B^  = 22.2/B^ 
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Equations  (13)  to  (20)  define  the  strumming  frequency  and  response 
amplitude  of  the  Seacon  delta  cables  to  a flow  field. 


VI.  DRAG  AMPLIFICATION 


In  Figure  7 (5,  13,  14,  15,  16,  17),  the  ratio  of  the  drag  coeffi 
cient  C on  a strumming  rigid  cylinder  to  the  drag  coefficient  on 
the  stationary  cylinder  is  plotted  as  a function  of  the  wake  response 
parameter  w defined  by 


Here,  Y is  the  amplitude  of  the  cylinder  vibration  and  tv  its  frequency 
The  solid  line  in  this  figure  is  given  by  the  formula 


which  is  a least  squares  fit  to  the  experimental  data  points  including 
the  constraint  C /C  = 1 when  w = 1. 


Experiments  by  Ramberg  and  Griffin  (14,  18)  have  shown  that  the 
near  wake  pattern  behind  a point  on  a strumming  flexible  cylinder  is 
basically  identical  to  the  near  wake  pattern  behind  a rigid  cylinder 
strumming  with  the  same  amplitude  and  frequency  as  the  point.  That  is, 
the  near  wake  pattern  behind  a point  on  a strumming  flexible  cylinder 
is  essentially  independent  of  the  spanwise  variation  in  amplitude  about 
that  point.  Using  this  information,  it  can  be  surmised  that  the  local 
drag  coefficient  Cp(s)  on  a strumming  flexible  cylinder  is  given  by 
equations  (22)  as 


where  the  local  wake  response  parameter  w (s)  is  defined  by  equation 
(21)  as  r 

w (s)  = [1  + 2Y(s)/d]  (in/ ii'  ) (24) 


In  particular,  for  a Seacon  delta  cable  strumming  in  its  n mode,  we 
have 

w (s)  - [1  + 2Y  (s)  /d]  (,„/«  ) (25) 


where  Y (s)/d  is  given  by  equation  (15) 


Finally,  we  obtain  ||ie  virtual  drag  coefficient  CpV  for  a delta 
cable  strumming  in  its  n mode  as 


CDV/CDO  = (l/L)  J [CD(S)/CD0]  ds  (26) 

It  is  this  virtual  steady  state  drag  coefficient  which  should  be  used 
in  array  motion  models. 

VII.  NUMERICAL  RESULTS 

In  Table  2,  the  calculated  values  of  C /C  for  each  of  the  three 
Seacon  II  delta  cables  are  listed  for  various  current  magnitudes  and 
angles  supplied  by  Albertsen  (9).  The  value  of  the  virtual  steady  drag 
coefficient  C due  to  strumming  is  frequently  1507.  to  2307.  larger  than 
the  value  of  the  nominal  stationary  drag  coefficient  C__.  These  large 
increases  in  the  steady  drag  coefficient  would  be  expected  to  have  a 
significant  effect  on  the  magnitude  of  the  array  motions  predicted  by 
the  array  motion  models  (2,3,4). 


In  Figure  8,  supplied  by  Albertsen  (19),  the  measured  motions  of 
a point  near  the  intersection  of  cables  1 and  3 on  the  delta  are  com- 
pared to  the  predicted  motions  during  a semi-diurnal  tidal  cycle.  The 
predicted  motions  were  obtained  using  the  Fortran  IV  program  DESADE 
(2) . The  magnitude  and  direction  (predominantly  from  the  southeast) 
of  the  current  during  the  cycle  can  be  inferred  from  the  displacement 
of  the  point  of  measurement  from  its  no  current  location  designated  by 
"X". 


The  measured  motion  of  the  point  is  shown  in  Figure  8 by  the  solid 
line.  The  predicted  motion  using  a constant  drag  coefficient  C = 

C = 1.55  is  given  by  the  dashed  line  with  circles  (--0--).  Finally, 
the  predicted  motion  using  C = (C  /C  ) C , with  C /C  calculated 

from  equation  (26)  and  = 1.55,  is  snownby  the  dashed  line  ( ). 

Excellent  agreement  between  the  predicted  and  measured  motions  is 
obtained  when  the  virtual  steady  drag  coefficients  are  employed.  It 
is  also  apparent,  for  the  larger  current  magnitudes,  that  the  displace- 
ments predicted  using  the  nominal  stationary  drag  coefficient  (C  = 

C = 1.55)  are  considerably  smaller  than  the  measured  displacements. 


VIII.  CONCLUSIONS 


The  most  recent  cable-strum  modeling  techniques  have  been  used  to 
calculate  the  vibration-induced  virtual  increases  in  the  steady  state 
drag  coefficients  of  the  Seacon  II  delta  cables  over  their  nominal 
stationary  values.  These  amplified  drag  coefficients  have  been  used 
in  the  array  motion  model  DESADE  (2)  to  predict  the  motions  of  the 
Seacon  II  cable  arrav.  Agreement  between  the  predicted  and  measured 
motions  is  excellent. 

The  results  of  this  work  clearly  show  the  importance  of  consider- 
ing strumming  and  its  resultant  drag  amplification  in  the  design  of 
mooring  arrays.  It  must  be  noted,  however,  that  much  research  remains 
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to  be  done.  At  the  present  time,  the  techniques  used  in  this  paper  for 
calculating  drag  amplification  are  limited  to  cables  in  a uniform  flow. 
Thus,  for  example,  nothing  can  yet  be  said  with  accuracy  about  any 
possible  drag  amplification  for  the  Seacon  II  anchor  legs  which j along 
their  lengths,  are  influenced  by  currents  that  vary  both  in  magnitude 
and  direction. 
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■I 

TABLE  1.  CHARACTERISTICS  OF  THE 

/ -C4-\ 

SEACON  II  DELTA 

m*  ( lb  s2  'l 

CABLES 

M*  ( Ik 

1 

V ft2  ) 

Mi  V 

Cable  1, 

T = 611 

lb. 

d 

= 0.712  in. 

1 

1.5 

0.0215 

3.56 

2 

999.5 

0.0215 

2.33 

3 

1001.0 

0.0215 

— 

Cable  2, 

T = 662 

lb. 

d 

= 0.715  in. 

1 

1.5 

0.0216 

3.56 

2 

333.0 

0.0216 

10.81 

3 

666.0 

0.0216 

10.83 

4 

999.5 

0.0216 

3.56 

5 

1001.0 

0.0216 

— 

Cable  3, 

T = 701 

lb, 

d 

= 0.625  in. 

1 

1.5 

0.0243 

2.33 

2 

76.5 

0.0243 

1.28 

3 

78.0 

0.0243 

0.00 

4 

81.5 

0.0152 

2.00 

5 

381.0 

0.0152 

1.44 

6 

556.0 

0.0152 

3.47 

7 

681.0 

0.0152 

1.44 

8 

776.0 

-0.0152 

3.34 

9 

856.0 

0.0152 

1.38 

10 

921.0 

0.0152 

2.28 

11 

923.0 

0.0152 

0.00 

12 

924.5 

0.0211 

1.28 

13 

999.5 

0.0211 

2.33 

14 

1001.0 

0.0211 

— 

ft  / 


Mass  plus  added  mass  of  water. 
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TABLE  2.  VIRTUAL  STEADY  STATE  DRAG  COEFFICIENTS 


V(ft/s) 

6v(deg) 

Cable  1 

Cable  2 

Cable  3 

.180 

162 

1.74 

1.78 

1.64 

.176 

159 

1.57 

1.30 

1.77 

.164 

154 

2.06 

1.89 

1.72 

.176 

142 

1.93 

1.80 

1.61 

.185 

152 

1.77 

1.71 

1.62 

.184 

143 

1.90 

1.73 

1.93 

.190 

141 

1.80 

1.68 

1.53 

.200 

134 

2.20 

1.65 

1.61 

.248 

120 

2.28 

1.51 

1.00 

.221 

115 

2.18 

1.73 

1.49 

.261 

123 

2.25 

1.96 

1.00 

.277 

141 

2.16 

1.71 

1.63 

.233 

102 

2.18 

1.89 

1.96 

.240 

102 

2.14 

1.91 

1.96 

.315 

98 

2.33 

1.64 

1.58 

.297 

96 

2.27 

1.79 

1.55 

.335 

84 

2.35 

1.27 

1.69 

.320 

84 

2.29 

1.32 

1.50 

.303 

78 

2.26 

1.70 

1.66 

.289 

76 

2.20 

1.77 

1.68 

.317 

68 

2.28 

1.18 

1.31 

.335 

70 

2.31 

1.00 

1.27 

.303 

62 

2.28 

1.00 

1.26 

.287 

51 

2.22 

1.22 

1.63 

.244 

58 

2.14 

1.00 

1.58 

.226 

41 

2.15 

1.00 

1.56 

.185 

18 

1.74 

1.36 

1.58 

.151 

31 

1.92 

1.00 

1.28 

.135 

117 

2.17 

1.48 

1.00 

CABLE  2 


.CABLE  3 


CABLE  1 


Fig.  5 — Relative  orientation  of  the  Seacon  II  delta  referred  to  an  x-y  base 
system.  The  current  at  the  delta  is  given  by  V and  the  angle  of  the  current 
with  respect  to  the  base  system  by  0V. 


'he  maximum  scaled  modal  response  amplitude  Amax  which  can  be  indui 
shedding  versus  the  response  parameter  SG  . The  solid  line  is  a least  squa 
eg end  for  the  data  points  can  be  found  in  Ref.  (7). 


Fig.  7 — The  ratio  of  the  drag  coefficient  CD  on  a strumming  rigid  circular  cylinder 
to  the  drag  coefficient  CDO  on  the  stationary  cylinder  as  a function  of  the  wake  re- 
sponse parameter  wr.  The  solid  line  is  a least  squares  fit  to  the  data  points.  Legend 
for  the  data  points:  o Griffin,  Skop,  and  Koopmann  (5);  + Griffin  and  Ramberg  (13); 
o Ramberg  and  Griffin  (14);  ® Tanida,  Okajima,  and  Watanabe  (15);  X Diana  and  Falco 
(16);  ♦ Meyers  (17). 


